On lower bounds of the solutions of some simple 
reaction-diffusion equations 
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Abstract 

The mean-field reaction-difi'usion equations of the diffusive pair-annihilation and triplett-annihi- 
lation processes are considered. A direct lower bound on the time-dependent mean particle-density 
is derived. The results are applied to the mean-field theory of the diffusive pair-contact process. 
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1 Introduction 



The interplay of reaction kinetics with particle diffusion has since a long time been a topic of intensive 
research. Here we shall consider the following reaction-diffusion equations 

dta{t,r) = Aa{t,r) - \a{t,rf (1.1) 

and 

dtb{t,r) = Abit,r)- libit, r f (1.2) 

where t G M+ and r G M" are time and space coordinates, A is the Laplacian and A, /i > are 
constant reaction rates. It is well-known that eq. (Il.lj) provides a mean-field description of the pair- 
annihilation process A + A ^ together with single-particle diffusion and similarly, p.2|) describes 
triplett annihilation A + A + A — > (we have rescaled the diffusion-constant to one). Still, the derivation 
of the long-time behaviour of solutions of such non-linear partial differential equations is not completely 
trivial. It is convenient to consider eq. ()1.1|1 inside a spatial domain f2 C M" and to define the mean 
densities 

a = a{t) = — dra{t,r) ; b = b{t) = — dr b{t,r) (1.3) 
l"l Jn l^^l Jn 

where \Q\ denotes the volume of Q. 

The most simplistic treatment of p.lj) simply suppresses the diffusion term, which leads to dtcis = 
—Xal and the solution as{t) = ao/{l + aoAt). Then, it may be asked to what extent this drastic 
simplification may be justified. This, and more generally the long-time behaviour of the space-time- 
dependent particle density a{t, r), has received a lot of attention. For example, a priori estimates such 
as the strong maximum principle may be invoked to obtain bounds V-{t) < a{t, r) < f+(t) where v±{t) 
satisfy the simplistic equation dfVj- = —Xv^ with the initial conditions ^^-(0) < a(0,r) < f+(0) ^21 
p. 94]. The latter conditions might however be difficult to meet for very 'rough' initial data a(0,r). 
For domains with a finite (and 'small') volume one may define an invariant region S and define 
a := rj — M, where rj is the principal eigenvalue of —A on and M = max^ iVAa^l. If a > 0, then 
it can be shown that a{t,r) converges exponentially fast with a characteristic time 1/a towards the 
solution as{t) \12\ p. 223]. However, the implied exponential approach need not hold any longer in 
spatially infinite regions. Then methods based on a scaling ansatz of the form a(t, r) = t~"/^/(rt~^/^) 
permit to extract the long-time asymptotics of the solution from phenomenological scaling [Sj or by 
using rigorous renormalization-group arguments This kind of argument can also be extended to 
systems of reaction-diffusion equations and in particular the reactions fronts in the two-component 
system kA + kB with initial conditions such that there is a reaction front between A-rich and B- 
rich regions. In n = 1 dimension, it was shown rigorously for k = 1 ^D] and /c > 4 j2] that the particle 
densities and the reaction front satisfy a multiscaling behaviour and that the convergence towards the 
scaling solutions is controlled by algebraically (and not exponentially) small corrections in t. On the 
diffusive scale, when \r\/\^ ^ 1, the problem essentially reduces to dta = Aa — |:2J. Reaction- 
diffusion systems of the form dta = Aa — ua — f{a) +g with a|j=o = <^(o)) <^\dn = 0, u > and suitable 
assumptions on / are reviewed in ^] . The long-time behaviour of the unique solution (roughly, on the 
Sobolev space W'^''^{Q) where g G L'^{Q) with q > max(2, n/2)) is described in terms of an attractor 
whose complexity can be analysed through its Kolmogorov entropy in great detail. 

Here we wish to present a simple direct estimate on the mean densities a{t) and b{t). We have (see 
section 2 for notations) 
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Theorem: (i) Let a G C^(]R+; W'^'^{Q)) be an (almost) everywhere non-negative solution of il.l]) . Let 
in addition be Va = on the boundary dQ. Then there is a positive and \Q\-independent constant X' 
such that the the mean density a satisfies 

- X'a{tf < dta{t) < -Aa(t)2 (1.4) 

(a) IfbE C^{W^; W'^'^{Q)) is an (almost) everywhere non-negative solution of / li.jjj) such that V6|^j^ = 
0, there is a positive and \ fl\-independent constant /i' such that 

- li'bitf < dtbit) < -iJibitf (1.5) 

For n < 3, the condition b E C^(]R+; W'^''^{n) n L'\n)) is sufficient. 

The Sobolev space Vr^'^(i7) may be too restrictive for unbounded domains. In section 2 we define a 
generalized space W'^'P{Q) for which the theorem stiU holds and which includes spatially homogeneous 
initial states. 

The upper bound in ()1.4j) has been known since a long time In particular, ()1.4pi.5|l imply for 
t sufficiently large (i.e. aoX't > 1 and 2bQfi't > 1, respectively, where oq > and bo > are the initial 
mean densities) the bounds 

^<ta(t)<i; <t^/^b(t) < ^= (1.6) 



It is admissible to take the infinite- volume limit oo. Eq. ()1.6|) is in full agreement with the results 

established earlier by different means [21 El El El 1121 • -^^^ similar upper bounds on ||a||p, (with 
1 < p < oo) which explicitly depend on the initial data, see [ZI- We point out that our derivation 
makes no explicit reference to the initial conditions (beyond the requirement a(t, r) > 0, b{t, r) > 0) 
and neither a scaling ansatz is needed. The bounds ()l.(i|l reproduce the expected mean-field scaling 
a{t) ~ t~^. For the mathematically oriented reader we recall that in low dimensions n < 2, the 
description of the diffusive pair-annihilation processes through a more microscopic approach such as a 
master equation (where fluctuations are taken into account) leads to a different long-time behaviour 
amicro(i) ~ t""/^, which has also been observed experimentally for n = 1, see jH| for a recent review. 
This manifests once again the character of equations such as (11.111 as mean-fleld approximations. 

The approach of a{t,r) towards the mean density can be described as follows. 

Corollary: Under the same conditions as in the theorem, there is a constant K' < 1 such that for 
times satisfying the conditions used in eq. il.f^) 

l^^dr {a{t,r)-a{t)f <^-t-' 

^J^dr {b{t,r)-b{tf<^.t-' (1.7) 



As an application, we consider the pair-contact process 2A — > 0, 2A 3A with single-particle 
diffusion (pcpd). We consider a domain Q C with the boundary condition Va|aQ = 0. The 
mean-fleld reaction-diffusion equation is 

dta{t, r) = DAa{t, r) + Aa(t, rf - /ia(t, rf (1.8) 
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with constants A G M and /i > 0. If the diffusion constant D = 0, a = a{t) evolves for A > towards a 
steady-state density a^o = X/fi, while aoo = for A < 0. It is known that 01 

r O (e-*/") ; if A > 
a{t) - Ooo ~ < t-^^^ ; if A = (1.9) 

[ ; if A < 

as t — cxD and where r > is a known constant. On the other hand, for a non- vanishing diffusion 
constant, we scale to D = 1 and have for the mean density 

Aa(t)2 - fi'a{tY < dta{t) < X'a{tf - fM{tf (1.10) 

for A > and 

- |A'| a{tf - fi'a{tf < dta{t) < -|A| a{tf - fM{tf (1.11) 

for A < 0, respectively and we can now let \Q\ oo, if we so desire. For A > 0, there is an active 
steady-state with density A//i' < ctoo < X'/fi but if A < 0, one has Ooo = 0. Furthermore, eq. (|1.9|) 
can be taken over. This proves the existence of a continuous steady-state transition at Ac = of the 
mean-field equation ()1.8|) . Finally, the approach of a{t,r) towards the mean density is according to 

^ r O (e~2*/-') ; if A > 

— / dr (a(t, r) - a{t)f <{ ; if A = (1.12) 

l^^l ''^ {t~^ ; if A < 

In n = 1 dimension, fluctuation effects create a very rich behaviour of the pcpd which is under active 
investigation, see jU] for a review. 

In section 2, we recall some inequalities which are needed in the proofs. The upper bounds in 
()1.4pi.5j) are derived in section 3 and in section 4, the lower bounds are obtained. 



2 Mathematical background 

We recall here some standard notations and some inequalities which will be needed in establishing the 
lower bound in (|1.4|1.5|) . We shall work with the p-norms, for 1 < p < oo 

dr\u{r)\p] (2.1) 



Denote by Lp{VL) the space of (equivalence classes of) functions with \\u\\p finite. Here and in the 
following VL C M". Furthermore, if VL has a boundary, it is assumed to be sufficiently smooth throughout. 
The space L°°(r2) is defined with respect to the supremum norm ||u||oo = ess supQ|u(r)|. For derivatives, 
we use the multiindex notation a = (ai, . . . , a„) G Nq and where |a| := ai + . . . + an- Then derivatives 
are denoted by 

d^u = -— — ^ (2.2) 

These derivatives can be taken to be weak (distributional) derivatives. The Sobolev space is 

W'''P{Q) := {u G LP{Q), d^u G LP(^]) for all |a| < k} (2.3) 

with its norm = X]|a|<fc 1 1'^"''^! Ip- Finally, C^{^y,W^'^{VL)) is the space of functions which are 

continuously differentiable with respect to t for all times < t < oo and whose values at any given t 
are in W''^P{yL). 
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After these preparations, we can state some known results which we need later. The first one is the 
Gagliardo-Nirenberg inequality, see [HI HI] • 



Lemma 1. For functions u G W^'^{VL)r\L'^{VL) with 1 < q < oc and for any integer < j < k, there is 
a constant C > such that 

\\d'u\\r<C\\u\\l-'\\d'u\\l (2.4) 

where 

--^- = + 0( 2.5 

r n q \p n J 

and if 1 < p < n / {k — j) , then j/k < 6 < 1. On the other hand, if 1 < p = n / {k — j) , only j/k < 6 < 1 
is allowed. 

We shall need two special cases of this. For n>2, wesetp = r = 2,/c = q'=l. Then j = and 
9 = n/{2 + n) G [|, !)• We get (here and in the following, we suppress the integration variable and 
write = J^dr) 

f r \4/(2+n) / . N n/(2+n) 

< Ci / |n| / \du\^] (2.6) 

Eq. ()2.6|) does not hold if n < 2. For 1 < n < 2, we set r = 2, p = g = A; = 1. Then j = 0, 

= n/2e [i, 1] and 

|2 



n 



u\' <C2{ \u\\ / \du\ (2.7) 



where the constants Ci^2 equal from Lemma 1. It can be checked from dimensional analysis that Ci 
and 6*2 are independent of 

Next, we quote an inequality due to Nirenberg, see 0. 

Lemma 2. For u G W'^''P{Vt), p > 1, there exists a constant Eq = eo{p,^) such that for any e with 
< £ < £o; there is a positive constant c = c{p, fl) such that 

[ dr \du\P <- [dr \uf + e / (2.8) 

\a\=2 

Finally, we quote Poincare's inequality, see fT]. Let Br{0) be the ball of radius R around the origin. 

Lemma 3. For any u G iy^'P(5H(0)) with 1 < p < oo, there is a positive constant Cp^ such that 

[ dr \u-u\P < C^f^|5^(0)|P/" [ dr \du\P (2.9) 

Jbr{0) JBr{Q) 

and the mean value u is defined in analogy with M.^) . 

It is sometimes desirable to consider spaces which are less restrictive than the spaces U'{VL). We 
define WiVL) as the space of (equivalence classes of) functions such that rripiu) := is finite. 

For unbounded domains (e.g. Vt = M") a limit procedure must be used in the definition of mp{u). We 
also set ^ 

Vr'^'P(fi) := |m G Lf{n), d"u G L^in) for all |a| < A;| (2.10) 

As an example, consider the function / : M — > M, a; hh^ f[x) = fo ^ 0. While / G Lp(M), since 
^pif) = |/o|^; clearly / ^ Lp{M.). Lemmas 2 and 3 readily extend to the space W'''P{fl). 
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3 The upper bound 



We briefly recall the proof of the upper bound in ()1.4|) . following pp. The mean density satisfies 

dta = / dr Aa — / dr 

\n\ \n\ 

If X p 

do- ■ Va- — / dr (3.1) 



1^1 Jdn 1^, „ 

where cr is a normal vector to the boundary dQ. The first term describes the flux of particles through 
the boundary and vanishes either in the limit of large volumes \Q\ ^ cxd or else if the boundary condition 
ValoQ = is imposed. Then dta = — Aa^ < — Aa^ by the Cauchy-Schwarz inequality. 

A similar argument works for triplett annihilation. By Holder's inequality, dtb = —fib^ < —fib^. 



4 The lower bound 

In order to obtain the lower bound in ()1.4|) . we recall from section 3 that dta = —\o?. The right-hand 
side is now estimated through the Gagliardo-Nirenberg inequality. We have to distinguish the cases 
n>2 and n < 2 and obtain from eqs. ()2.6|) and ()2.7p 



I -AC,N-(/.a)'/^^-^"^(/jVan"/(^^") ;if.>2 



(4.1) 



where for n < 2 the Cauchy-Schwarz inequality was used again. Next, we need an upper estimate for 
|Vap, which is provided by the following 

Proposition: For u G l^^'^(i7) there is a constant c > and an e^, such that < < oo and that for 
all e < one has 

< — [ \u\P (4.2) 



n ^ Jn 

Proof: If ^|q|=2 l*^"^!^ ~ proposition holds true trivially, because of ()2.8|) . We can thus 

suppose that J2\a\=2 In > 0. Next, the function f{x) := A/x + Bx, where A,B are positive 

constants, has an absolute minimum at x=k = ^A/B. For X < x^,, the first term dominates over the 
second and f{x) < 2 A/x for all x < x*. We apply this to the inequality ()2.8|) of Lemma 2. The 
right-hand side is minimal if e = e^,, where 



p 



1/2 



with the £0 of Lemma 2. Then the assertion follows. q.e.d. 
The extension to Vr^'^(f2) is immediate. 

Therefore, setting p = 2 and appealing to dimensional analysis, there is a positive constant K > 
such that 
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Jn 
Jn 



< 2K\n\'^-^/''a^ (4.3) 

From the eqs. ()4.1|1 it follows 

C^a > I _^c2i2Ky^/^ a2 ; if n < 2 ^^'^^ 
This is exactly the form asserted in the theorem and we can identify the effective reaction rate 

r ACi(2irr/(2+") ;ifn>2 
^ \ AC2(2ir)"/2 ; if n<2 ^^-^^ 

Remark: the last bound in ()4.3|) might be improved by restricting to a ball around the origin and 
applying the Poincare inequality. If C^^^K < 1/2, one gets jVap < K/{1 - KC'^^) l^^l^-^/^a^. For 
n = 1, the bounds might be further sharpened with the help of the inequality 



f \du\P <ce~^^^^ ( f \u\] +e [ \d\\P 
Jn \Jn J Jn 



(4.6) 



where n{p) = — (p — 3 + 1 /p) /p and which can be proven for p > 1 through a slight generalization of the 
proof [H] of the inequality ()2.8p of Lemma 2. Since /i(p) > if 1 < p < = (3 + -\/5)/2, the method of 
the proposition can for example be used in the p = 1 and p = 2 cases. 

The lower bound in ()1.5|) is proved similarly. First, we set A; = g = 1 and r = 3 in Lemma 1. Then 
j = 0. For 1 < n < 3/2, we set p = 1, find ^ = 2n/3 G [2/3, 1] and 

Next, for 3/2 < n < 3, we set p = 3/2, find 6 = 2n/{n + 3) G [2/3, 1] and 

/ n \ 3{3-n)/(n+3) / ^ n 4n/(n+3) 

{ \u\] ( / \du\''/^) (4.8) 



'n \Jn / \Jn / 

Finally, for n > 3, we set p = 3, find = 2n/(2n + 3) G [2/3, 1) and 

„ / f \ 9/(2n.+3) / ^ N 2n/(2n+3) 

/ \u\'<tJ \u\) / \dun (4.9) 

Jn \Jn J \Jn J 

Dimensional analysis shows that Fi 2,3 (which stand for of Lemma 1) are independent of For 
the further analysis of ()4.9|) . the proposition above states that there is a positive constant Kj, such that 
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|V6p 



n 



6^ 



+36 / |6-6| 



1-3 



(4.10) 



Now, from dtb = —fib^, we obtain, using first eqns. ()4.7|4.8p4.9|) . then the Cauchy-Schwarz inequal- 
ity and the following consequence of Haider's inequality ^ |V6|=^/2 < |^]|i/4 (J^ | v6|2)=^/^ and finally 
eqns. (|4.3I4.1()|1 



-f^^iM-' (/J V6|) 



, 2n 



> 



-fir,\Q\-' ^J^^)3(3-n)/(3+n) ^.^^ | 3/^^^ 4n/(3+n) 

-/ir2|firi+"/("+3) (/^6)=^(=^-")/(=^+") (^|V6|2)'"/^'+"^ 



if 1 < n < 3/2 
if 3/2 < n < 3 
if 3 < n 

if 1 < n < 3/2 
if 3/2 < n < 3 
if 3 < n 



1-3 



> 



-/iri(2ir)"6 
-/ir2(2ir)3"/("+3)6^ 



,1-3 



if 1 < n < 3/2 
if 3/2 < n < 3 
if 3 < 



(4.11) 



-/ir3(5ir3)'"/^'"+^^6' 

which is the form asserted in p.5|l . We identify the effective reaction rate 



/iri(2ir)" 
/ir2(2i^)3"/("+3) 

|2n/(2n+3) 



if 1 < n < 3/2 
if 3/2 < n < 3 
if 3 < n 



(4.12) 



This completes the proof of the theorem. q.e.d. 

The bounds (jl.6j) are established by direct integration, since dta(t) ^ —Xa{ty implies a{t) ^ ao/(l + 
aoXt) from which the assertion is immediate, under the stated condition on t and for gq > 0. The bounds 
for b{t) follow similarly. 



We now prove the corollary. First, \Q\ ^ J^{a — a)^ < \Q\ 



-1 r 772 _ 772 



a 2 < (At) Next, we try to 



improve this bound by letting = Bji{0) be a ball around the origin and apply Poincare's inequahty 

^ / (a-a)2<42)|^|-i+2/n f \^a\'<2KC^'^a' 



(2) 

and therefore K' = min(l, 2KCp ). The estimates for b are obtained similarly. 



q.e.d. 



Remark: for the PCPD with D = 1, eqns. p.lOll.ll|) are obtained by applying the results of the 
theorem separately to the two terms on the right-hand side of eqn. p.8|) . Eq. ()1.12|) is obtained by 
replacing a 1— *■ a — a^o in the corollary and then using p.9|) . Let ati^) '■= a{t,r). For sufficiently long 
times, at G W^'^{W), but if A > 0, a* ^ W^'^{W). 
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